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ABS TRACT 
Solitary  waves  in  running  polytropic  gases  of  finite 
or  infinite  height  are  studied.  Explicit  expressions  for 
the  critical  speed  and  the  solution  of  the  solitary  waves 
are  obtained  by  a  perturbation  scheme  applied  to  the  non- 
linear equations .   It  is  found  that  internal  waves  may 
also  be  observed  in  a  running  gas,  and  change  of  wave  tj^e 
may  occur  in  the  medium.  The  results  obtained  are  expressed 
in  terms  of  simple  integrals  of  density  and  velocity 
distributions  of  the  medium  in  the  state  of  equilibrium, 
and  conform  with  solutions  for  liquids  of  constant 
density  as  special  cases. 
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1.  Introduction. 

Since  1834  when  Scott  Russel  made  his  first  observation 
of  a  solitary  \mve   on  a  stream,  much  theoretical  vjork  has 
been  done  on  the  problem  of  solitary  waves  in  a  liquid  of 
constant  density.  We  refer  the  reader  to  Stoker's  book, 
"Water  Waves"  [1],  for  a  detailed  discussion  of  the  problem 
and  the  references  cited  there.   In  recent  years,  there  has 
been  growing  interest  in  the  study  of  solitary  waves  in  a 
stratified  medium.  An  interesting  account  of  the  recent 
developments  of  the  solitary  wave  problem  has  been  given 
in  [2]  (Peters  and  Stoker,  i960).  A  variant  of  Friedrichs ' 
procedure  is  used  by  Peters  and  Stoker  to  solve  the  problem 
of  solitary  waves  in  liquids  having  non-constant  density, 
and  the  so-called  "internal  waves"  are  found  on  the  inter- 
face of  a  two-layer  liquid.  Recently,  Benjamin  also  solved 
the  problem  of  solitary  waves  on  a  stream  with  an  arbitrary 
distribution  of  vorticity  [5].   In  all  the  previous  work, 
the  medium  considered  is  always  assumed  to  be  incompressible, 
and  the  height  of  the  medium  finite.  For  the  work  done  here, 
we  shall  be  concerned  with  the  problem  of  solitary  leaves  in 
a  running  polytropic  gas  of  finite  or  infinite  height.  The 
main  difficulty  encountered  in  our  problem  is  that  the 
density  as  well  as  the  pressure  vanishes  on  the  free  surface, 
which  may  be  at  infinity,  and  the  equations  become  singular 
there.   In  order  to  overcome  this  difficulty,  v;e  use  a 
transformation  established  by  Peters  and  Stoker  to  express 
all  the  functions  in  terms  of  the  horizontal  distance  x  and 
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the  stream  function  ip   as  nevj  Independent  variables,  and 

the  limiting  values  of  all  the  functions  as  the  Image  of 

the  free  boundary  In  the  x,V''~plane  Is  approached  from 

below  are  then  defined  as  the  values  of  the  function  there. 

Furthermore,  since  ^  Is  a  function  of  vertical  distance  t, 

In  the  state  of  equilibrium,  we  may  eventually  formulate 

the  problem  in  the  x,C-plane  as  Benjamin  first  did  for  the 

Incompressible  case.  However,  the  approach  applied  here  is 

different  from  that  in  [ 3] . 

We  formulate  our  problem  in  §2.   The  solitary  wave 

solution  is  obtained  in  §5.  The  coefficients  of  the 

solitary  wave  equation  are  reduced  to  simple  Integrals 

of  the  density  distribution  p^  and  the  velocity  distribution 

u  in  the  state  of  equilibrium.   It  is  shown  in  §4  that  the 
00 

known  results  for  liquids  of  constant  density  are  special 
cases  of  our  solution.   Furthermore,  the  so-called  "neutral 
stream  lines"  are  found  and  some  criteria  for  the  existence 
of  such  lines  are  established.  As  a  consequence,  the 
internal  waves  may  also  be  observed  in  a  running  gas.  The 
appendices  contain  the  details  of  derivations  which  we  have 
omitted  in  the  main  text.  The  results  obtained  here  are  only 
for  the  case  of  a  polytropic  gas  with  arbitrary  equilibrium 
velocity  profile,  but  it  is  feasible  to  extend  the  method 
to  gases  with  more  general  pressure  and  density  relations . 
We  hope  that  the  research  reported  here,  if  relevant,  may 
be  of  some  use  to  the  study  of  certain  geophysical  phenomena 
In  the  upper  atmosphere. 
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This  report  is  intended  to  supersede  Part  I  and  Part  III 
of  an  earlier  report  [4],  the  results  of  which  are  only- 
special  cases  of  the  solution  given  below.  At  the  same  time, 
we  also  take  this  opportunity  to  correct  a  few  errors  commit- 
ted in  the  computation  of  the  coefficients  of  the  solitary 
wave  equation  given  in  [4]. 

Solitary  Wave  in  a  Running  Gas 
2.  Formulation  of  the  Problem. 

We  consider  a  layer  of  a  polytropic  compressible  medium 
supported  by  a  rigid  plane  bottom  and  having  a  free  surface 
which  may  or  may  not  be  at  infinity.  The  pressure  at  the 
free  surface  is  assumed  to  be  zero  and  there  are  no  geometric 
constraints.   In  the  state  of  equilibrium  the  velocity  profile 
as  a  function  of  the  vertical  distance  only  is  assumed  to  be 
given.  A  cross  section  of  the  medium  is  a  horizontal  strip 
in  the  upper  half  plane  (Fig.  1).  It  is  supposed  that  a 
wave  of  permanent  type  moving  to  the  left  with  constant 
velocity  c  has  been  generated  in  the  medium  as  a  result  of 
some  disturbance.  We  choose  a  coordinate  system  moving  with 
the  wave  such  that  the  x-axis  coincides  with  the  bottom  and 
the  y-axis  passes  through  the  crest  of  a  wave  of  elevation 
or  the  trough  of  a  wave  of  depression  and  is  positive 
upward  (Fig.  1).  With  respect  to  the  coordinate  system  the 
wave  is  stationary  and  the  medium  at  infinity  is  in  the  state 
of  equilibrium. 

The  equations  governing  a  tvjo-dimensional  steady  flow 
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of  a  polytroplc  compressible  medium  are 


dx      ay 


•«  Su  ,  ~  c)U      1  Sp 

^  -S-x  ■*■  ^  ^  =  -  ^  ^ 
~  Sv  .  ~  3v     ^   1  Sp 

^  ^  +  ^  -^  =  -  e  -  ^  Bf 

Pq     Po 

where  p(x,y)  andp{x,y)  are  respectively  the  pressure  and 
density,  u(x,y),  v(x,y)  the  horizontal  and  the  vertical 
velocity  component,  g  Is  the  gravitational  constant,  and 
Pq,  Pq  are  the  reference  pressure  and  density.   In  the 
state  of  equilibrium,  we  have 

~  n      n-1  Cil/n-l 


(2.2) 


e 


for  n  >   1,   0   <  C   <  Cs=  ~i  h. 


Po 

for  n=l,   0<^<C=co 


r\j  i\, 


where  h  =  p^/gpQ,     ^  Is  the  vertical  distance  at  equilibrium 
and       ,  pQ  is  the  value  of  p^  at  c,   =  0 . 

Denote  by  tl/{x,y)    the  stream  function  such  that 

^y    =    PU     ,  ^j^    =    -    PV     . 

The  mass  flux  across  any  vertical  plane  from  y  =  0  to  the 
free  surface  y  =  y„  per  unit  breadth  is  given  by 
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'/'(x^Yg)  -  -^(x,0)  =  J     ^^  u^  dy  =  7g  , 

0 
where  we  may  set  ^(x,0)  =  0. 

Let 

^(x,y)  =  7 

where  0  <  y  <  y„,  0  <  7  <  7^,  and  assume  that  for  any  value 
of  7  in  0  <  7  <  7  there  exists  a  unique  solution  of  the 
above  equation  for  y  such  that 

y  =  f(x,7)  . 

Here  the  bar  notation  is  used  to  denote  a  quantity  as  a 
function  of  x  and  7.   If  we  choose  f,  p,  p,  and  u  as  dependent 
variables,  then  we  have  from  (2.1),  for  0  _<  7  <  7„,  -00  <  x  <  «>, 

^x  =  ^V^x  -*-  V7 

^^xx  +  ^x^x  =  -  PS^7  "  ^7 
(2.3) 

pUf^  =1,     V  =  uf^ 

I-   =  (i-)"  ,     n  >  1, 

Pq   Po 

subject  to  the  boundary  conditions 

f(x,0)  =  0,     p(x,7s)  =  0.^^^ 

(•^'It  is  noted  here  that  the  transformation  from  x,y-plane 
to  x,7-plane  breaks  down  as  7=73  since  fy=co   there.  Hence 
we  define  the  limit  of  a  functionas  7  ->  73  as  the  value 
of  the  function  at  7  =  7  .   In  fact,  7  =  7g  is  a  branch 
line  in  the  x,7-plane.  For  a  general  discussion  of  the 
mappings,  c.f.  [5]. 
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However,  In  order  to  make  u  and  p_  appear  explicitly  in 

^  CO  •  00 

the  above  equations,  it  will  be  more  convenient  to  use 
the  vertical  distance  C  at  equilibrium  as  an  independent 
variable.  The  relation  between  y   and  C  Is  given  by 

y  =/pco  ^00  ^y  ' 

0 
and  assuming  that  67/ 6£,  >  0,    for  0  <  ^  <  ^g,  we  have 

si  _  1  H 
00 
where 

and  the  hat  notation  denotes  a  function  of  x,^.   Since  p^(C) 
is  always  positive  for  0  <  ^  <  h^  where  h^  =  f(x,7g),  we 
suppose  that 


1 

C3 


U  (O  7^   0 


for  0  <  C  <  Cg- 


By  using  x,C  as  independent  variables  (2.3)  becomes, 


for  0  <  C  <  C^  -00  <  x  <  00, 
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^<«^^^xx  +  ^x^x^  =  -  PS^C  -  PC 

^   g      CO  '  X 

I-  =  (&-)''  ,         n  >  1  , 

Po   Po 

subject  to  the  boundary  conditions 

?(x,0)  =  0  ,  p{x,^g)  =  0  . 

Now  we  introduce  the  following  dimensionless  variables 


^ 

= 

X 

h   ' 

-n 

= 

^ 

f 

= 

f 
h     ' 

p^ 

Pco 

-  -^     > 

0 

V 

u 

V    =    —    , 

u  =  —  , 

c    ' 

c    ' 

A. 

/\ 

p 
P  ~  'v^  2   * 

p  =  ^  , 

Po° 

Po    A 

^=^   , 

G     =        ~   , 

c 

P  »c 

^«, 

^s 

"c«    =    c        ' 

'Is   -   h      ' 

A*      <V 


where  h  =  Pq/ePq*   and   (2.4)   become,    for  0   <  n   <  o    ,    -co  <  |    <  co. 


G  U/,   =  -   f  Pfi    +  f^p 

00    I  ',-j^^  ^^T) 


G    (uffi>    +  u,,  f.  )    =   -   pT^f^    -  p,. 


puf^  =  G     ,  V  =  uf^    , 

p  =  Ap"   ,  n  >   1   , 


subject  to  the  boundary  conditions 


f(C,0)  =  0  ,  P(e,r,g)  =  0 
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3"   Nonlinear  Theory.   Solitary  Wave  Solution. 

It  is  assumed  that  a  solitary  wave  moves  with  a  speed 
such  that  ?v  =  gh/c   is  near  some  positive  value  S, ,    which  will 
be  determined  later.   The  equations  (2.6)  can  he   written  as, 
for  0<r]<ri,   -oo<|<oo, 

G  u^  =  -  f  p^  +  f.p 


(3.1) 


G^(uf^,+u^f^)  =  (,(;-A)pf^-.gpf^  -  p.^^ 

puf^  =  G^  ,         V  =  uf^  , 
p  =  -(i-A)p^  +  -2p"  ,   n  >  1  . 


We  let 


s  =  i  -  A 
and  introduce  a  new  variable 

Then  in  terms  of  -q,  a,  (3.1)  become, 

for  0   <  -q    <  -q    ,-w   <   a   <  «>, 
s 


(3.2) 


G   u^  =   -    f  p^   +  f  D 

00  a  q^cr         at] 

eG    (uf     +u  f   )   =  epf     -  ipf     - 

puf^   =  G  ,                     V  =  euf _  , 

p  =   -£p  +iJp      ,               n>l. 


subject  to  the  boundary  conditions 

f(cj,0)  =  0  ,     p(a,qg)  =  0. 
Suppose  that  all  the  dependent  variables  can  be  expanded 
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in integral  powers  of  e,  i.e. 

(3.5)  ^{cs,r\,e)    =   >   s  (k  (a,n)  , 

k=0    ^ 

where  ^   stands  for  p,  p,  u  and  f.  Substitution  of  (50) 

in  (5*2)  will  j^ield  a  sequence  of  equations  and  boundary 

conditions  for  the  successive  approximations  by  letting  the 

coefficients  of  like  powers  of  s  be  equal. 

The  equations  for  the  zero-th  order  approximation  are, 

for  0  <  ri  <  iT  ,  -M  <  a  <  CO, 
s 


(3.4) 


.u.f...  =  G^ 


Pq  =  Jlp^   ,        n  >  1, 

subject  to 

fQ(a,0)  =   0  ,     p(cT,Tig)  =  0. 

(5'4)  is  a  set  of  non-linear  equations,  the  solution  of 
which  cannot  be  found  by  any  standard  method.  However, 
we  may  assume  Pq  =  p^  since  the  solitary  wave  solution  is 
essentially  a  small  perturbation  over  the  equilibrium  state 
at  a  =  -CO.  Based  upon  this  ansatz,  we  find  the  solution 
for  the  zero-th  approximation  as  follows. 
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(3.5) 


Po  =  P«.  ' 


f 0  =  ^. 


0       00 


The  equations  for  the  first  order  approximation  are, 

0<ri<ri  ,  -oo<a<oo, 
s 

^«>^la  =  -  %Pla  +  ^laPori  ' 

0  =  Po^op  -  -^Po^m  +  Pi^oii^  -Pl^  ' 

Pi^o^^'on  ■"  Po^i^oti  ■*•  Po^o^iTi  =  ^  > 


(5.6) 


subject  to 


Pi  =  -^np^PQ   -  Pq 


f^(a,0)  =  0  ,  Pi(a,r|g)  =  0  , 


where  use  is  made  of  the  fact  that  pQ,pQ,fQ  and  Uq  are 

functions  of  •,-;  only.  Elimination  of  f,,p,,  and  u,  from  (;^.6) 

and  making  use  of  (5.5)  yield  the  equation 

I 

(5-7)  (ufpi.„).,-  =  Y   '"'Pina'  • 

I  00 

The  derivation  of  (5.7)  is  deferred  to  Appendix  I.  From 

(3.7)  we  have 

(5.8)  Pl^^  =  --^  ^a'(a) 

CO 

and  since  Pi-(ct,t|  )  =  0,  we  also  obtain 


,X   <  n 


t         q    -   ,-.0 


''■   "   0"^ 


rX 


ci  ~- 


-■    =    rQ 


^^   -    (.i- 


•  o  *> '? 


[a 


.   ( 


n  r. 


>•) 


(V^^) 


liK^'-a;'?      ,1  vj  ,  0-? 


wj. 


-jfila  SriB 
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(5.9)  P-Lp  =  ^a'(a)F(n)  , 

where  a '(a)  Is  an  arbitrary  f-unctlon  of  c   and 

Prom  the  last  equation  of  (^^.6)  we  have 

1  ^-n+1 
Pla  ~  Jin   PQ   Pla 

=  -k   Po''"''^a'(a)P(ri)  , 
and 

Pi  =^  p-"'^^a{a)P(ri)  +b(ri) 

where  ^^;e  assiime  a  (a)  ^  0.  Since  o  tends  to  p^  as  a  —>   -<» 
and  we  also  choose  a(-t»)  =  0,  the  above  equation  shovjs  that 
b(Ti)  =0  and 

Pi  =  i  Po"''^a(a)F(-a). 


It  follows  that 


Pi  =  '"PiPo   -  Po 

=  ^a(a)P(ri)  -  ^'^Pq 


Prom  (3.'^')>  (5.6),  (3.8)  and  (5-9),  we  find  that 

^la  =  Po^'^^^^"^^%T.a-^Vla^ 
=  -  a'(a)[-J3'^+p"^P(ri)]  , 


and 


f 3_  =  -  a(a)[-e'-^+  p;^^P(r|)]  +b^(n) 


brtiS 


S''BCl 


)   10 


Ivr. 


ri.  - 


T>I 


,      ^}-{o) 


■Kf 


Cj  oi 


9fi     ,0    .-.v     (03-).. 
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c''''S    -    ff:)'^v(o)::?.   - 


i.e^U  brirl  sw    .(>:)    b;::::,   C:  .1)    J-.r) 


or         r 

■oL  '••.•';       r::i  '■•  ■      '  0  ^       dI 
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Slnce  f,  — >  0  as  0  ->  -co,  \-ie   have  b^(r, )  =  0  and 

f ^  =  -  a(a)[-^""^  +  p;^F(r|)]  . 
By  applying  the  condition  f^(cT,0)  =0,  we  see  that 

(5.10)  .e  =  p^-^Co)  . 

In  summary,  the  solutions  for  the  first  order  approximation 
are 

p^  =  iJa(a)F(n)  -  ^"■'■Pq  , 

Pi  =^  Po"''^a(a)F(r^)  , 

(3.11) 

f^  =  -  a(a)[-.e"^+p;^F(Ti)]  , 

u.  =  -  -(^^ 


1  "    ^«>   ' 

00 

To  determine  the  function  a(a)  we  must  proceed  to  the 
equations  for  the  second  order  approximation,  vjhich  are,  for 
0  <  T)  <  -n  ,  -w  <  a  <  cr, 

^J^'Jlaa^    =    ^P0^1ri-^l%P^--'^^P0^2r,+Pl%+P2S^-P2^    ' 
(3.12) 

P0^2S  +  PO^CS  ^  P2^0^0i.  -^  Po^l^lTi  +  Pl^oS  "^  Pl'-l-^On   =   °    ' 

»   n-1     „   n-1  ,  n(n-l)  „  2^n-2 
subject  to  the  conditions 


-ci- 


b:\.0    0   3=    {  ;■;)  |-' 


fl 


90  2   "^v:    ,0  ^ 


(j)^'"'^^  =■■   I 


r_ 


t     .^f:''"3.    -     (  :",j'^(t))ji:..    =     r<5 


fo)^ 


^-  -  =  I'^ 


.0   >   T) 


.C  4- 


~>-. 


q)        ~      ^ 


«    r:oq"(.o^c--^.r'irr 


q)    -    '■ 


?     sr     -   •'■       "r     ff    n  +        "t     (!     '^  •'-        '■"     f?    n  4-        't    ff 


r  r 


"o'l'''^  -  Q':^'on«'   ^^  cjli 


sriol  J.L5:'tOO    9- id'    c- 


'^'  ^=  („r./-')r.cj 


-    (0,d)„ 


-ih. 


Prom  the  solutions  obtained  for  the  zero-th  and  first 

approximations,  vje  can  obtain  the  following  equations 

2 

u  Po  ^      1 

(—7 )-,  =  —  g4(a,Ti)  , 


(5.15) 


where 


Pco 


^20=    ^P0Ti)'^(^'"^"^P2-aa-^P2a^0T  +  ^^"^S2^-ei^   ' 


gl(a,r))  =  -  p^^g^^^  +  f^^p^^  , 
g2(cT,il)  =  G^u^f^^^  +  G^u^  +  &p^ul   , 
gjCa^n)  =  -  ufp3_f3_^  +  p^u^  , 

■*■  Po-aPo^^-^C^g2a-*-^^^~Sa^ 
+  P0TiTiP0-n(S2a--Si)  +  P^Z^.^    • 

The  derivation  of  (3.I5)  x\;ill  be  given  In  Appendix  I. 

Now  using  the  condition  ^oA^}'^)   =  0,  we  obtain  from  (5. 10) 

and  (3.15)  that 

(5.1^)   P(0)u^(0)p2^^((7,0)  +  P2^(cT,0)  +  P(0)g2^(a,0)  -  g^(a,0)  =  0 

In  order  to  eliminate  Pp   (a,0)  and  Pp  (cr,0)  from  {'^.Ih),   we 

must  make  use  of  the  first  equation  of  (5.I5).   Multiplying 

both  sides  of  the  equation  by  F(t))  and  then  Integrating  with 

respect  to  -o  from  0  to  ri  ,  v;e  have,  by  Integration  by  parts, 

s 

2  ^s 

ujDp        -n     p  gi,(CT,i^. ') 


_i;r-. 


crprix"    briB 


^riol:)£';po 


LLC'l   o 


evr    , 


t    '. 'T-.r^j 


*■         :        CO    ■ 


-  =   { 


(01. C)   mo'T'r  nlBacfo  ev:    .0  -   ^■•■ .  o )  ^ -.1  noljibnoo  o- 
■  ■  ■     t>'.^ 


C   =    {.-.o)r-    -    (0 


t:)^Qa(''')^  +    (O.'o)-^c;'j  -•'    f  ■■..'^). 


3r-i;vIc:i::MMil      .f"X.C)    "o  HOi::fj:r;/'j.    •■3i.f:V    sr'-t  ':.o  e\ 


LI    - 


(•n)1^- 


'.;) 


-15- 


It  can  be  shovin  that  P('Oufo"-^Po„^  ->  0  as  ri  ->  v^   If 
u^(r]    )   j^  0,   and  since  Pg^^^^'^'s^   ^  ^*    ^^^  above  equation 


becomes 


S 


-  P(0)u^(0)p2      (a,0)    -  P2^(<J,0)    =    /^    F(ti  '  )g4(a,P  •  )p'"^(ii  ' )    dr^ 


;ari 

0 


It  follows  from  (3-1^)  that 
^s 


-  y'  P(il')Sii(a,n')p;^(T^,')  dn'  +  F(0)g2^(a,0)  -  s^{a,0)  =  0 
0 

This  gives  the  condition  to  determ.ine  the  function  a(o-). 

By  some  lengthy  but  straightforward  calculations,  we  finally 

obtain  the  equation 


tit 


where 


mQa  ( a )  +  m^a ' ( a ) a ( a )  +  mga ' ( a )  =  0 


•p 
s 

0         ^1  00   ■   •  CO   '  '00 

b 

^•s 
m^  =  -  5  /^  .gp^(-p)u;^^(Ti)  dn.  -  p^(0)P(0)  , 


(3.16)  0 

mg  =  P(0)  , 

^  =  P"-^(0)  . 

The  derivation  of  these  results  will  be  given  in  Appendix  II, 
Suppose  that  none  of  m^,  m,  and  mg  are  equal  to  zero,  and  we 
also  impose  the  conditions 

a'(-oo)  =  a"(-co)  =  a'(0)  =  0  , 
Then  the  solution  for  a(a)  is 


!~! 


rfb  (• 


{■o)b  I'C  ii:c;:' :    -■ '-}   snirrrcs-'tv^   or^  ncirji: 


fn-'B^:-  -h    ( 


•'I 


d) 


iti 


^"O^    )    ^:--     )''i'^ 


.-'■•T 
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r 


(0)^" 


.i  a 


i       -0 


:'•)•<?     ^     (r>-.-)"S 
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AssujTie  that  the  successive  approximations  up  to  the 

first  order  give  a  sufficiently  accurate  representation  of 
a  solitary  vjave.  We  have,  in  terms  of  the  independent 
variables  x,  and  t,,    for  0<^<ri^,   -co<x<  -f<o. 


P  =  Pn?"-(^-A)PpF^.)  -^  sech^  ^V-  iT  ^--^) 


^0'^co-^*-'^^t'0''^';^  ^  '^""^'   2hy   m^ 


f  ^  hri  -  (.C-A)h[,r^-p;^P(-n)]  -j^  sech^  ^  y-  ^  (.C-A)  , 
(5.17)  u  ^  cu^  +  (^-?.)c(u^)  1  —^  sech^  |j-  /-  -^  (.e-A)   , 


"CO  •'••"•  OS'     m,       2h  y   mQ 


v 


s  (,5-^)'/2  jou„[ri  - 


CO    .      m^ 


If  the  equation  of  state  for  a  perfect  gas 

p  =  RpT 

is  used,  where  R  is  the  gas  constant,  and  T  the  absolute 
temperature,  then 


A. 


Rpo         ''Po  1     I-   0 


,'.'--r    > 
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f  r.,  ■■'  ", 


,-.!P      ■ 
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4.  Discussion. 

Prom  (3.1")  we  see  that  itIq  is  always  negative;  if 
jg~"^-p^-^((-il)  s  0,  then  itIq  =  0,  and  also  f^  h  0  as  seen  from 

(3.11),  and  no  solitary  wave  exists.   It  is  evident  that 

n'r\         Poo^'H) 

mp  =  P(0)   =    /  — ^ drj   is  always  positive.     As  observed 

^0    uii-T])  nip 

from  the  solution  for  a  (a), =■  (i-A)  must  be  positive  to 

m^ 


ensure  the  existence  of  a  solitary  wave,  and  it  follows  that 
&'7\  >   0  and  c^  >  ^''''gh  =  P(0)  ~  .  In  Appendix.  II,  it  will 

be  shown  that  m,  is  always  negative.  The  wave  type  of  the 

-1  -1 
solitary  vjave  is  determined  by  the  sign  of  &     -0^  F(-o)  as 

seen  from  the  expression  for  f  in  (3. 17).  We  have  a  wave 

—1  —1 
of  elevation  (depression)  If  <l      -p^  F(n)  >  0  (<  0).   Let 

us  write  i'^^-o  7{t])   as  — [P(0)o  -F(-n )  ]  .   It  is  easily  seen 

•  00    '        0  'CO 

that  g(-n)  =  P{0)p^-F(ri)  =  0  when  ■(]   =  0,-f\^.      In  order  to 

determine  whether  there  are  any  zeros  of  g{r))    in  0  <  -n  '^  "^.^t 

we  established  some  criteria  by  considering  the  poliits  of 

intersection  of  p  F{0)  and  P(n)  other  than  r\  =   0,'n_.  We 

shall  first  find  the  slopes  of  p^F(O)  and  F{t\)   at  'o  =  0,  r]=j\^^ 

Now  p^('n)F(0)  =  -^  p'J'-^^i'.^FiO),    F^  {^)    =  -  p^u^^. 

At  n  =  0,  ,  n 

p^(0)P{0)  =  -^  P(0)  , 


P'(0)  =  -  u;^(0)  , 


>  n  >  1  ; 


and  as  ri  — >  r   , 
s 


p^('n)F(0)  ->  -00  ,  00  >  n  >  2  , 

->  -  I  F(0)  ,  n  =  2  , 

— >  0  ,  2  >  n  >  1  , 

P  •  (ri )  ->  0  ,  n  >  1  . 


1-:^    .Jl    t> 
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1      '      <     i 


-      ^ 
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_      -       :  I 


,1      <      S 
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Therefore,  for  n  >  1,  P'('q)  <  p^P(O)  near  ti  =  "„•   If 
P'(0)  >  p^(0)F(0),  then  P(r,)  >  p^P(O)  near  t^  =  0.   It 
follows  that  there  must  exist  at  least  one  zero  of  q{'.]) 
in  0  <  ri  <  -n  .  Vfe  may  write  the  condition  P'(0)  >  p'(0)P(0) 
for  n  >  1  as 

o 

u  P(0)  >  n  ,         for  n  >  1. 

For  n  =  1,  and  -n  near  ri^  =  w 


If  we  assume  that 


P^(t^.)P(0)  =  -  p^(ri)P(0)  , 


F'(r:)  ^  -  oJr^U'Jir..^)       . 


U^('0g)P(O)    >    1      {<   1), 


in  addition  to  u^(0)F(0)  >  1  (<  1),  there  also  exists  at 
least  one  zero  of  P(0)p^-P('p)  in  0  <  r,  <  n  .  The  stream 
lines  in  0  <  ■■•■   <  ^„   along  wl  ich  f,  =  0  are  called  "neutral 

A 

Stream  lines."  Since  f-,  =  0  at  o  =  0,  there  must  exist  one 

extremum  of  i^    in  0  <  n  <  v     if  a  neutral  line  occurs. 
1  3 

Therefore^  "internal  waves"  may  also  be  observed  in  a  running 
gas.  Since 


for  n  >  1, 


for  n  =  1, 


lljii  p;^P(ri)   =  0    ; 


lim  p;^F(v-)   =  ufi-.)    , 
lim  [olh{r^,)]^  =  0    , 
lim  [p"-'"P(o)]   ^   =   2  u"^(n    )u'(n„)    , 

11->TlS  •   ' 


n 
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f,  is  an  absolute  maximum  at  n  =11  ^or   n  >  1,  and  is  a 
relative  maxiraum  (minimum)  at  p  =  r;  if  u'(ii  )  ^  0  and 

S        00   t> 

We  shall  compute  a  simple  example,  pedagogical  as  it 
may  be,  to  illustrate  what  we  have  discussed.   Let  n  =  2,  and 

by  (2.2).  IJe  assume  that 

u^  =  (l+^-o)"^/^  . 


00 


Since 


^;^(^'^  rpoo^c^  dii  =  3  >  2 , 


0 

there  exists  a  neutral  stream  line  in  0  <  r]  <  2.  Mow 

1    "^   r;   P 

which  has  zeros  at  -n  =  0,  1 '2  and  2.   The  minimum  of 
p(0)-p^-'-F(r;)  equal  to  -I/I6  occurs  at  tj  =  1/4,  and  the 
maximum  equal  to  5  occurs  at  the  free  surface  ii  ==  2.  The 
solitary  wave  is  a  wave  of  depression  in  0  <  ri  <  1/2,  and 
becomes  one  of  elevation  in  1/2  <  ti  _<  2. 

In  the  following  we  consider  some  special  cases  of  our 
solution. 
Case  I.   p  =  constant. 

'  CO 

This  case,  in  fact,  corresponds  to  the  one  of  n  — >  <». 
We  have  p  =  1  and  n  =  1  by  normalization.  Then  from  (3.1(i), 


.- ; 


.0   = 


:;-  o^ 


,  CJ     <- 
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0  1 

0  1^ 

0 


As  shov?n  in  Appendix  II, 
1  rs   -n" 
niQ  =  P.  J  f  ful^{^)   ulby)    u;2(.r,M)    dT)   dr.'    d^i" 
0   0   0 

-    r  '^^!('0    u"^(n')    d-o    diV    . 
0  0 
We  find  from  (3.17)    that,   at  \^,  =  1, 

1 

f(x,i)  =  h  +  h(?v~^-.r^)(Ar^)[  fui^  d  ]"^ 

0 


ch^^  [A(?.-^-.r^)/(-mo)]^/2 


se..     2h 

Since  'KZ"'^  =  1  +  0(£),    the  above  equation  becomes 

1  1 

f(x,l)   =  h+h(A'^-iJ"^)[?v   fnl^  dTil"^  sech-|^  [(A''^-,r^)/(-r\)r, 

0 

which  is  the  result  obtained  by  Benjamin  for  the  solitary  wave 

on  running  water  if  we  identify  Tv  with  g  in  the  equations  (32) 

and  {y4)   of  [J]. 

Case  II.  u  =  constant. 

CO 

without  loss  of  generality,  we  take  u^  =  1.  Making  use 
of  (2.2),  we  have  from  {3'!''^) 


'\D 


c~> 


r„ 


1 


0 


C: 


,   ^'2.  -  gin 


-i; 


!!  '.s    I  ~r- 


5    '  ifc    fr&    ( 


1~<- 


;5      ( 


'"•     i 


^I    -   '-.   ::5    ,"Bi:i 


J. 


•    ](■'   U)!' 


--.).  .. 


U 


rn~)\c'^-^--^-^)^]  ^ 


r 
c 


?,sniO':-ac  c;C'i1"Bupe  ©■/oi'b  g;M    .(3)0-!-  x 

I 


l-).'^"-^"/^)]   ;%  :;o£)3   -^"[fTo 


si;?  nx 


;■   9vr 


'U 


e-A£.:i   e\!    ^-^,- 


(■!.-)    wnl 


-21- 


-2(n-l)^   ^^^ 
"b  ~  (2n-lH3n-2)  ' 


IHt  = 


-3n+l 


1  "    n       ' 

nip  =  1  ,    ^  =  1  . 

As  n  — >  03,  ni-,  — >  -1/3,  m,  — >  -3.  ¥e  have  the  same  solitary 
wave  equation  for  a  liquid  of  constant  density  obtained  by 

Peters  and  Stoker  if  we  set  r  =  0,  i  =  1,  and  p=  1  in  the 

equation  (^1.20)  of  [21.  As  n  =  1,  mQ=  0,  and  no  solitary 
vjave  exists. 


^■'■^The  coefficients  m^,  m^  and  m^   obtained  below  should 
replace  those  given  in  Part  I  of  [4]. 
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Appendix  I. 
The  equations  for  the  first  approximation  are 

(1.2)  0  =  pof^^    -  .C(pof;i_^   -i-  p-^fo^)    -  Pi^^      , 

(1.3)  PiU^fQ     +  pQU^fQ     +  p^UQf^     =  0      , 


/-r   I,  \  ,'  n-1         n 

(I.^)  Px  =  -^■^Pi?o       '  PQ      ' 


subject  to 

fx(a,0)   =  0   ,  Px(a,ns)   =  0    . 

From  (3.5)  and  (1.3).^  we  have 

where  f^^  =  1.  Multiplying  the  left-hand  side  of  (1.2)  by 


u  ,  and  substitulng  the  above  equation  for  (Pi+Poo^in  ^ '  ^^ 
obtain 


(1.6)  0  =  G<«  +  -G^^l  -  ^fPln  • 

Then  by  differentiating  (1.6)  with  respect  to  a  and  ti 
successively,  it  is  found  that 
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Maklng  use  of  (1.7),  we  obtain  from  (I.l)  and  (I. 3)  that 


Novj  differentiating  both  sides  of  (I.l)  with  respect  to  r\ 
and  making  use  of  (1.9)  and  (I. 10),  we  have 

^"^f^^Plna^  =  -  Pla^i  -^  Po-nP5^(-^'^Plc^-  Pia^ 

+  Por^.^Po^-'^^^lcrri+Pla^  ' 
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Appendix  II. 

We  recall  that  the  condition  for  a (a)  is 
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Here  we  note  in  passing  that  all  the  integrals  are 

proper  since  F(rj)  =  O(p^)  as  ri  — >  rig. 

The  coefficients  itIq,  m^  and  m^   can  be  greatly  simplified 

by  showing  that  many  terms  are,  in  fact,  equal  but  of 
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By  integration  by  parts  of  the  second  member  on  the 
right-hand  side  of  the  above  equation,  we  have 
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In  the  follovjing,  we  shall  show  that  the  sum  of  the  last 
four  integrals  on  the  right-hand  side  of  the  above  expression 
is  equal  to  zero.  First  v.'e  note  that 
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by  integration  by  parts.      In  the  follov,?jLng,   we  shall  show 
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Director 

Grad.  Sch.  of  Aeronautical  Engrg. 

Cornell  University 

Ithaca,  New  York  (1) 


Chief,  Bureau  of  Aeronautics 
Department  of  the  Navy 
Washington  25,  D.  C. 
Att:   Research  Division       (1) 

Chief,  Bureau  of  Ordnance 
Department  of  the  Navy 
Washington  25,  D.  C. 
Att:   Research  and  Dev.  Div.   (1) 

Office  of  Ordnance  Research 
Department  of  the  Arroy 
Washington  25,  D.  C.         (1) 

Commander 

Air  Res.  and  Development  Command 

P.  0.  Box  1395 

Baltimore  I8,  Maryland 

Att:   Fluid  Mechanics  Div.    (1) 

Director  of  Research 

National  Advisory  Comm,  for  Aeron. 

1724  F  Street,  N.  W. 

VJashlngton  25,  D.  C.         (l) 

Director 

Langley  Aeronautical  Laboratory 
National  Advisory  Comm.  for  Aeron. 
Langley  Field,  Virginia       (1) 

Director 

National  Bureau  of  Standards 

Washington  25,  D.  C. 

Att:   Fluid  Mechanics  Section  (1) 

Professor  Richard  Courant 

Courant  Inst,  of  Math.  Sciences 

New  York  University 

251  Mercer  Street 

New  York,  New  York  10012     (1) 

Professor  G.  Kuertl 
Dept.  of  Mechanical  Engineering 
Case  Institute  of  Technology 
Cleveland,  Ohio  (1) 
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Chief,  Bureau  of  Ships 
Department  of  the  Navy- 
Washington  25,  D.  C, 
Att:   Research  Division       (1) 
Code  420  Prelim.  Design  (1) 


Commander 

Naval  Ordnance  Test  Station 
3202  E.  Foothill  Blvd. 
Pasadena,  California 


(1) 


Commanding  Officer  and  Director 

David  Taylor  Model  Basin 

Washington  7*  D.  C. 

Att:  Hydromechanics  Lab.  (l) 
Hydrodynamics  Div.  (1) 
Library  (1) 

(1) 


Ship  Division 


California  Inst,  of  Technology 
Hydrodynamics  Laboratory 
Pasadena  4,  California        (1) 

Professor  A.  T.  Ippen 
Hydrodynamics  Laboratory 
Massachusetts  Inst,  of  Technology 
Cambridge  59,  Massachusetts   (1) 

Dr.  Hunter  Rouse,  Director 

Iowa  Inst,  of  Hydraulic  Research 

State  University  of  Iowa 

Iowa  City,  Iowa  (1) 

Stevens  Institute  of  Technology 

Experimental  Towing  Tank 

711  Hudson  Street 

Hoboken,  New  Jersey  (1) 

Dr.  L.  G.  Straub 
St.  Anthony  Falls  Hydraulic  Lab. 
University  of  Minnesota 
Minneapolis  14,  Minnesota     (1) 

Dr.  G.  H.  Hickox 
Engineering  Experiment  Station 
University  of  Tennessee 
Khowvllle,  Tennessee         (1) 


Chief  of  Naval  Research 
Department  of  the  Navy 
Washington  25,  D.  C. 
Att:   Code  4l6 
Code  460 
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Chief,  Bureau  of  Yards  and  Docks 
Department  of  the  Navy 
Washington  25,  D.  C, 
Att:   Research  Division 


Hydrographer 
Department  of  the  Navy 
Washington  25,  D.  C. 

Director 

Waterways  Experiment  Station 

Box  651 

Vicksburg,  Mississippi 


(1) 


(1) 


(1) 


Office  of  the  Chief  of  Engineers 
Department  of  the  Army 
Gravelly  Point 
Washington  25,  D.  C.  (1) 

Beach  Erosion  Board 

U.  S.  Army  Corps  of  Engineers 

V/ashlngton  25,    D.    C.  (1) 


Commissioner 

Bureau  of  Reclamation 

Washington  25,    D.    C. 


(1) 


Dr.  G.  H.  Keulegan 
National  Hydraulic  Laboratory 
National  Bureau  of  Standards 
Washington  25,  D.  C.         (1) 

Brown  University 

Graduate  Div.  of  Applied  Math. 

Providence  12,  Rhode  Island   (1) 

California  Inst,  of  Technology 
Hydrodynamics  Laboratory 
Pasadena  4,  California 
Att:   Professor  M.  S.  Plesset  (1) 
Professor  V.  A.  Vanonl   (1) 
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Mr.  C.  A,  Gongwer 

Aerojet  General  Corporation 

6552  N.  Irwindale  Avenue 

Asusa,  California  (1) 

Professor  M.  L.  Albertson 
Department  of  Civil  Engineering 
Colorado  A.  and  M.  College 
Fort  Collins,  Colorado        (1) 

Professor  G.  Birkhoff 
Department  of  Mathematics 
Harvard  University 
Cambridge  58,  Massachusetts   (1) 

Massachusetts  Inst,  of  Technology 
Department  of  Naval  Architecture 
Cambridge  59*  Massachusetts   (1) 

Dr.  R.  R.  Revelle 

Soripps  Inst,  of  Oceanography 

La  Jolla,  California         (1) 

Stanford  University 

Applied  Math,  and  Stat.  Laboratory 

Stanford,  California         (1) 

Professor  J.  W.  Johnson 
Fluid  Mechanics  Laboratory 
University  of  California 
Berkeley  4,  California       (1) 

Professor  H.  A.  Einstein 
Department  of  Engineering 
Universitv  of  California 
Berkeley  4,  California        (1) 

Dean  K.  L.  Schoenherr 

College  of  Engineering 

University  of  Notre  Dame 

Notre  Dame,  Indiana  (1) 

Director 

Woods  Hole  Oceanographic  Institute 

Woods  Hole,  Massachusetts     (1) 

Hydraulics  Laboratory 

Michigan  State  College 

East  Lansing,  Michigan 

Att:   Professor  R.  R.  Henry   (1) 
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